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First and second sound speeds are calculated for a uniform superfluid gas of fermi atoms as a 
function of temperature, density and interaction strength. The second sound speed is of particular 
interest as it is a clear signal of a superfluid component and it determines the critical temperature. 
The sound modes and their dependence on density, scattering length and temperature are calculated 
in the BCS, molecular BEC and unitarity limits and a smooth crossover is extrapolated. It is found 
that first and second sound undergo avoided crossing on the BEC side due to mixing. Consequently, 
they are detectable at crossover both as density and thermal waves in traps. 
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I. INTRODUCTION 

Recent experiments probe systems of strongly inter- 
acting Fermi and Bose atoms and molecular bosons at 
low temperatures. Evidence for a continuous crossover 
between BCS to a molecular BEC superfluid is found 
by measurements of expansion^ S S S 0> collective 
modes 00) RF spectroscopy [3, thermal energies and 
entropies [||, and correlations [l(| ■ The next generation 
of experiments will measure sound velocities in these sys- 
tems. Second sound has been studied in He 3 , He 4 and 
BCS superfluids (see, e.g., 0,0). Second sound is 
particularly interesting as it is a clear signal of a super- 
fluid component and determines the critical temperature 
|l3j. The dependence of the sound speeds on tempera- 
ture, density and interaction strength will constrain the 
equation of state and models. 

Eagles and Leggett described the smooth 
crossover at zero temperature by a mean field gap equa- 
tion. Extensions to finite temperature especially around 
the superfluid critical temperature can be found in 
E E El E 13, |H El E| . The various crossover mod- 
els differ in the inclusion of selfenergies, whether one or 
two channels are employed, etc. Recently, the crossover 
has been calculated in detail by quantum Monte Carlo 
(QMC) |25| at zero temperature. The experiments 
are generally compatible with most these calculations 
and confirm universal behavior |l.'il |26| in the unitar- 
ity limit as well as superfluidity at high temperatures 

m mm mm. 

The purpose of this work is to calculate sound ve- 
locities in the superfluid phase and demonstrate that 
first and second sound undergoes avoided crossing at the 
BEC-BCS crossover in uniform systems of strongly inter- 
acting Fermi atoms. Second sound reveals the presence 
of a superfluid and its critical temperature. The sound 
modes provide details about the equation of state at the 
BEC-BCS crossover and constrain crossover models. 

The two-body interaction range is assumed to be short 
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as compared to the s-wave scattering length a and the 
interparticle spacing or kp , where n = k F /3n 2 is the 
density of the Fermi atoms with two equally populated 
spin states. The physics can be expressed in terms the 
variable x = l/(akp), which varies from — oo in the di- 
lute BCS limit through the unitarity limit x = at the 
Feshbach resonance to x — > +oo in the molecular BEC 
limit. 

The sound modes may be measured in experiments 
as density waves in the expanding clouds. Density or 
thermal perturbations in the trapped clouds may be in- 
troduced in the center of the cloud, for example, by a 
pulsed laser before expansion, or by other means |29j| . 
It is important that the two sound modes mix as will 
shown below because both modes will be excited by ei- 
ther a density or thermal perturbation. As the sound 
modes travel through the cloud the density and temper- 
ature may vary. It may therefore be nontrivial to extract 
the sound speeds as function of density, temperature and 
scattering length, and useful to have a model calculation 
of both sound speeds. 

The manuscript is organized such that the sound 
modes are first described in general, and then calculated 
in the BEC, BCS and unitarity limits respectively. The 
three limits are then approximately matched assuming a 
smooth crossover. Finally, summary and conclusions are 
given. 

II. SOUND MODES 

We shall in the following assume that the thermal ex- 
citations collide frequently and therefore can be assumed 
in local thermal equilibrium. This is the collisional limit 
where hydrodynamics apply. The transition from zero to 
first sound as systems chan ge f rom collisionless to colli- 
sional is described in , e.g ., [30( and the hydrodynamics 
of superfluids in j^, yj- For a Fermi liquid the colli- 
sion time r oc T™ 2 can become large at low temperaq- 
tures. The hydrodynamic limit is therefore limited to 
long wavelength (or equivalently low frequency oj) sound 
modes such that tlo <C 1. This can always be achieved in 
a bulk system but not necessarily in a finite system where 
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collective frequencies u> are finite. However, it was found 
in 0, 0, that axial and radial collective modes can be 
described by hydrodynamics at temperatures well below 
the superfluid transition. Thus atoms in traps appear 
collisional in a wide region around the unitarity limit. 
Yet the collisionless limit is expected in traps at very low 
temperatures [34$ . 

The presence of a normal and superfluid components 
leads to two sound modes in the collisional limit referred 
to as first and second sound. Their velocities U\ and u 2 
are given by the positive and negative solutions respec- 
tively of mm 



(i) 



The thermodynamic quantities entering are the adiabatic 
c| = (dP/dp)s, the isothermal Cj, = (dP/dp)?, and 
the "thermal" c\ = p s s 2 T / p n cv sound speed squared. 
The latter also acts as a coupling or mixing term. The 
difference between the adiabatic and isothermal sound 
speed squared can also be expressed as 
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FIG. 1: First and second sound speeds for temperatures 
T/T c — 0.25, 0.5, 0.75. The isothermal sound speed ct is cal- 
culated from the Leggett crossover model. Sound speeds are 
calculated using entropies from Eq. Q on the BEC side, Eq. 
11211 on the BCS side and from the measured entropy in the 
unitarity limit at x = (see text). 
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Here, p = p n + p s = nm are the total, p n the normal and 
p s the superfluid mass densities, s = S/p the entropy per 
unit mass and cv — T(ds/dT) p the specific heat per unit 
mass. 

The isothermal sound speed at zero temperature is in 
both the hydrodynamic limit and for a superfluid gas 
given by 
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where vp — hkp/m, (3{x) — Ei nt /Ekm is the ratio of 
interaction to kinetic energy, j3' = df3/dx, etc. QMC 
calculations 0E1 find 0(x = 0) = -0.57 in the unitar- 
ity limit at zero temperature. The continuous crossover 
found in QMC and crossover models is confirmed exper- 
imentally for thepressure, chemical potential, (3, collec- 
tive modes 0, 0> @ an d recently also for the entropy . 

In the dilute BCS limit and at low temperature ct — 
Vf\J (1 + (2/Tr)akp)/3. In the unitarity limit ct = 
v F y/(l + /3(0))/3 ~ 0.37v F . In the dilute BEC limit 
ct = \J (tt /2)h 2 naM Jrn? , where (Zm is the molecular scat- 
tering length. According to QMC [2^| and four-body 
calculations j^ij «m — 0.62a, whereas in the Leggett and 
several other crossover models au = 2a. 

We are mainly interested in the superfluid state at tem- 
peratures below the critical temperature, which gener- 
ally is less than that in the BEC limit, T < T C <T® EC = 
0.218-Ep, the temperatures will be much smaller than the 
Fermi energy Ep (in units where ks = 1). At such low 



temperature the isothermal sound speed is given by Eq. 
(PI in both the BCS and BEC limits. In the crossover 
model of Leggett the isothermal sound speed is also to a 
good approximation given by the zero temperature value 
for T<T C . We shall therefore use Eq. @ for c T in the 
following. 

The entropy and mass densities have been calculated 
in the BEC [H H3 and BCS H3 limits to leading or- 
der in akp. One should, however, be careful in extrap- 
olating to the unitarity limit. For example a normal 
BEC is depleted as ajv/ — > +oo and eventually quenched 
|38| whereas in crossover models the condensate remains 
although it changes character from BEC to BCS su- 
perfluid. Likewise, the collective mode frequencies in 
trapped Fermi atoms increase for a normal BEC as the 
unitarity limit is approached [3^ . Crossover models [3^ 
and experiments jy, Q > however, find the opposite for a 
molecular BEC. The BEC and BCS limits will now be 
described followed up with a discussion on the extension 
and interpolation towards the unitarity limit in various 
crossover models. 

III. BEC LIMIT 

Crossover models and calculations find different molec- 
ular scattering length in the BEC limit as mentioned 
above. Also the effective mass M in the dispersion re- 
lation E q = q 2 /2M differs between the molecular BEC 
limit where M = 2m and the unitarity limit, where M 
becomes very large indicating that the dispersion rela- 
tion is not quadratic but rather linear as a Bogoliubov 
phonon. For example, the pair susceptibility Xj which 
enters selfenergics, is calculated with different combina- 
tions of bare (Go) and full (G) Green's functions. With 
X = GG [ItI |20| M increases monotonically as the uni- 
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modes caused by the mixing term in Eq. The mix- 
ing acts similar to level repulsion between otherwise de- 
generate states in quantum mechanical systems. Instead 
of crossing the two sound speeds avoid crossing by an 
amount proportional to the mixing term, and first sound 
turns continuously into second sound and visa versa. As 
the avoided crossing occurs predominantly on the BEC 
side it can be described by using the BEC mass densities 
and entropies in Eq. The resulting sound modes 

are shown in Fig. I with ct and T c from the Leggett 
crossover model which has au — 2a. For comparison the 
sound modes are calculated in Fig. 2 with ct from QMC 
PH . which has om = 0.62a. In both cases the entropy of 
Eq. with M = 2m is assumed. 



FIG. 2: As Fig.l but with c T and au = 0.62a from QMC. 
Furthermore, induced interactions are included on the BCS 
side which reduce T c . 



tarity limit is approached from the molecular BEC limit, 
whereas with \ — Go Go [la \lM M first decreases and 
then increases as x — ► 0+. With x — GqG ^3 M 
first increase, then decrease and finally increase again 
as a; — > 0+. The critical temperature has the opposite 
behavior as M because T c cx 1/M. 

A normal BEC has a linear Bogoliubov dispersion re- 
lation E q = qcT for temperatures below T* = 2naMn/M 
corresponding to i^(om /a)(T c /T). For oj/ = 0.62a and 
T ~ T c /2 the quadratic dispersion relation is thus lim- 
ited to x^jl. Similar results are found in several crossover 
model [l9l I23I ] , where the molecular gas can on the BEC 
side (x^jl) be described approximately by a dispersion 
relation E q — q 2 /2M, with an effective mass M ~ 2m. 

The corresponding entropy per particle is then given 
by the standard result 



S_ 

N 



3/2 



(4) 



where a = 5<(5/2)/3C(3/2) = 0.856 and the crit- 
ical temperature is T c = 2Trh 2 (n/2((3/2)) 2 / 3 /M = 
0.218Ep(2m/M). The normal mass density is on the 
BEC side p n = p(T /T c ) 3 / 2 . From the entropy and mass 
densities on the BEC side we obtain for the mixing term 
c 2 = a(T/M) Ps /p and c| - c 2 T = a(T/M)p n /p. 

In the BEC limit we thus obtain the classical result for 
the first sound speed 



u\ = \ a 



T 
M 



(5) 



and the standard second sound speed u\ = c^ps/p 
inh 2 aMPs/M 3 , or 



u 2 
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As we approach the unitarity limit we observe that 
avoided crossing takes place between the two sound 



IV. BCS LIMIT 

The thermodynamic quantities and sound speeds in 
BCS limit have been calculated in (3?J from the mean 
field BCS equations of Leggett 0] 
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with coupling strength g 
h 2 k 2 /2m - p, E k = y/e' 2 + A 2 and A(T) the gap. The 
thermal distribution function is = (exp(Ef./T) + 
With the equation for number density conservation 
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the gap and chemical potential can be calculated as func- 
tion of density, temperature and interaction strength. At 
zero temperature the last terms in Eqs. Q and ijHJ van- 
ish and the gap is 



A 



Ef exp(7r/2a/c_F) . 



(9) 



Gorkov |40| included induced interactions, which reduce 
the gap to 



A = (2/e) 7/3 ^exp( 7 r/2afc jF ) 



(10) 



Induced interaction can be included in the gap equation 



by replacing a 



- 2fc F ln(4e)/37r on the left 



is corrected but 
. In both cases 



hand side. Hereby, not only the gaj 
also the thermodynamic potential 3 
T c = (7/tt)A = 0.567A . 

The thermodynamic functions can be calculated from 
the thermodynamic potential per volume Q = —P. We 
make the standard assumption that the Hartree-Fock 
terms in the superfluid Q s and normal state f2 n ther- 
modynamic potentials are the same. The difference is 
then given in terms of the pairing coupling as 



Va'^ 
dA' 



(11) 
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The first order temperature correction to the ther- 
modynamic potential in the normal phase is f2„ = 
-A(0)tt 2 T 2 /3, where N(0) = mk F /2n 2 is the level den- 
sity. 

Inserting the coupling of Eq. Q into the thermody- 
namic potential l|llf> and again exploiting that A <C p, 
it reduces to 



-N(0) 



In j - AT J de k ln(l - f k 



(12) 



The crossover model thus arrives at the standard ex- 
pression for f2 s and therefore also the standard entropy 
density S s — —(dti s /dT)v,n, and specific heat C s = 
T{dS s / dT)v p, in superfluid phase. In the normal phase 
S n = C n = N(0)2ir 2 T/3. Near T c : S s /S n (T) = 1 - (1 + 
0(1 - T/T c ) and C s /C n (T c ) = £ - 3.77(1 - T/T c ). Here, 
£ = 1 + 12/7C(3) ~ 2.43 is the ratio of superfluid specific 
heat just below T c to the normal one just above T c . At 
low temperatures S s = N(0) v / 2ttA 3 ,/T cxp(— Ao/T) and 

C s = N(0)y/2nA*/T3exp(-&o/T). 

Finally, we need the superfluid (London) density 



n 1 



where n = n s + n n is the total density. At low tem- 
peratures n s /n = 1 — y 27r Ao /T exp(— Ao /T) whereas 
n s /n = 2(1 — T/T c ) near T c . 

We now have all the thermodynamic quantities avail- 
able for calculating the sound speeds in the BCS limit. 
Since thermal effects are small we find that u\ — cs — ct ■ 
The second sound is a pure thermal wave with velocity 



S 2 T 
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The dependence on temperature and x is shown in Fig. 
3 as well as in Figs. 1 and 2 for three different tempera- 
tures. 

At low temperatures the second sound is linear in tem- 
perature 
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Around T ~ 0.7T C the second sound speed has a broad 
maximum of 
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Near the critical temperature T c — T <C T c 
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The characteristic u 2 oc yl — T/T c behavior near T c is 
the same as in the BEC limit, Eq. ©, but with a differ- 
ent prefactor. 




FIG. 3: Thermodynamic quantities for a superfluid as calcu- 
lated in the BCS limit vs. T/T c . Shown are the gap A(T)/A , 
entropy S s (T)/S n (T), specific heat C s {T)/C n {T)£, and su- 
perfluid density n s (T)/n. The second sound is plotted in 
units of vfT c /Ef- 



V. UNITARITY LIMIT AND CROSSOVER 

At crossover the entropy receives contributions from 
thermal Fermi atoms and molecular bosons. The en- 
tropies and sound speeds reflect the underlying dispersion 
relation for low- lying excitations. On the BEC side for 
T^T* the quadratic dispersion relation E q = H 2 q 2 /2M 
leads to an entropy S oc T 3 ' 2 , whereas for T^ET* the lin- 
ear Bogoliubov phonon E q = hqcx yields S oc T 3 . On 
the BCS side the gap in the quasiparticle energy leads to 
a more complicated entropy as discussed above. 

Generally, for a power law dispersion relation E q = q a 
the entropy and specific heat scale as S oc Cy oc T 3 / Q 
and the normal mass density p n = p(T/T c ) 5 /" -1 . Con- 



-l 

2-2/c 



sequently cg-cf, oc T A l a+1 and c| oc (p s / p)T 2 - 2 / a . The 
resulting temperature dependence of the sound velocities 
follows from Eq. QJ. 

Model calculations [H, EE EE E3, EE 123, [13, E3i 

find a smooth crossover of most thermodynamic quanti- 
ties as function of density and scattering length at zero 
and finite temperature except at T c , where superfluidity 
vanishes. We will therefore extrapolate the sound modes 
from the BCS and BEC limits towards crossover such 
that they are continuous around crossover and meet the 
sound speeds in the unitarity limit, which now will be 
calculated independently. 

Experimentally, the entropy has recently been mea- 
sured and found to scale as S/N ~ 20(T/T F ) 153 in the 
unitarity limit at temperatures below T c ~ 0.27 Ep 0- 
In the unitarity limit the entropy is independent of scat- 
tering length due to universality EH as °t ner thermody- 
namic quantities j2(J. Consequently, when S/N oc T 3 / 2 , 
the entropy density is independent of density and is there- 
fore the same in bulk as in a trap at temperatures suffi- 
ciently low as compared to the Fermi temperature. Both 
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the measured T c and entropy are then slightly larger in 
the unitarity limit than for the pure BEC of Eq. 
With these adjustments of the parameters in Eq. J2J, 
the sound velocities are calculated in the unitarity limit 
as in the BEC limit above as shown in Figs. 1 and 2 
for x = 0. To guide the eye the unitarity limit point at 
x = is connected smoothly in the crossover regime to 
the BEC limit, x > 1, and the BCS limit, x < -1. 

The sound velocities in Figs. 1 and 2 are well behaved 
in the sense that both the BEC and BCS limits extrap- 
olate continuously towards the unitarity limit. On the 
BCS side first sound is basically given by the isothermal 
sound speed, and the second sound increases exponen- 
tially with x as T c = (j/tt)A as given by Eqs. © and 
(|10|l . Around the unitarity limit their dependence on 
temperature and x = 1/kpa is more complicated because 
the two sound modes mix and the temperature depen- 
dence of the entropy is different from the BEC limit. On 
the BEC side the curves again extrapolates continuously 
towards the BEC limit, however, with avoided crossing 
of the two sound modes. 

It should be noted that the curves in Figs. 1 and 2 are 
plotted for constant T/T c . On the BEC side T c is almost 
constant from the unitarity to the BEC limit, and the 
curves are therefore basically for constant temperature. 
On the BCS side, however, T c decrease exponentially as 
T c oc exp(7rai/2) (see Eqs.© and ltTU)l L and therefore 
the temperatures also decrease. If the temperature were 
held fixed, superfluidity and the second sound would thus 
vanishes below a certain x where T > T c . 

The specific heat is a continuous function of tempera- 
ture around T c only in the BEC limit. It is discontinu- 
ous in the BCS limit by a factor CJCn = 2.43. Recent 
calculations and experimental fits [l9l | at x = 0.11 find 
C s /C n — 2.51 ± 0.05 at T c . These experiments indicate 
that the discontinuity in the heat capacity has a maxi- 
mum crossing over from the BCS value towards zero in 
the BEC limit. T c also has a maximum at crossover in 
several crossover models 0, 0, ^| . 

The second sound speed squared generally vanishes 
with the superfluid density at the critical temperature, 



i.e. C2 oc y/T c — T. This is observed in the BCS limit of 
Eq. U3)), in the BEC limit of Eq. © and also in the 
unitarity limit |13|. 

The sound speeds are continuous as function of x at 
crossover as long as T/T c is fixed as shown in Figs. 1 
and 2. As function of temperature first sound depends 
only little on temperature whereas second sound van- 
ishes at T c and the specific heat is discontinuous. Similar 
temperature dependency were found in Refs. [2ll l4l| at 
crossover. 



VI. SUMMARY 

The first and second sounds have been calculated in 
the molecular BEC, the BCS and the unitarity limits. It 
has been argued that the sound modes interpolate con- 
tinuously between these limits, and that avoided crossing 
takes place on the BEC side due to mixing. 

The mixing has the consequence that both sound 
modes can be excited and detected both as density and 
thermal waves in traps. In contrast the second sound is 
a pure thermal wave in the BCS and BEC Bogoliubov 
limits. In finite traps the sound waves could be excited 
in the center and propagate towards the surface so that 
the local density decreases. Consequently, \x\ = l/\a\kp 
increases and the sound velocities change as shown in 
Figs. 1 and 2. 

The second sound mode reveals superfluidity and the 
critical temperature. Measurements of both sound modes 
constrain the equation of state, entropies, specific heat, 
etc. and thereby crossover models, which predict very dif- 
ferent behavior of, e.g., the molecular scattering length 
in the BEC limit and effective molecular mass, critical 
temperature, etc. at crossover. The crossover and mix- 
ing of the sound modes will also be possible to study in 
mixtures of bosons and fermions. A number of different 
superfluids may also be created in optical lattices 
(with corresponding Mott transitions) and possibly also 
supersolids. 



[1] K. M. O'Hara, S. L. Hemmer, M. E. Gehm, S. R. 

Granade, J. E. Thomas, Science 298 (2002) 2179. M. 

E. Gehm et al., Phys. Rev. A 68, 011401 (2003); 
lcond-m at/0304633]. 
[2] M. Greiner, C. A. Regal, D.S. Jin, Nature 426, 537 

(2003); Phys. Rev. Lett. 92, 040403 (2004). 
[3] T. Bourdel et al., Phys. Rev. Lett. 91, 020402 (2003); 
[4] M.W. Zwierlein et al, Phys. Rev. Lett. 91, 250401 

(2003); S. Gupta et al, Science 300, 47 (2003). 
[5] S. Jochim et al., Phys. Rev. Lett. 91, 240402 (2003); M. 

Bartenstein et al., ibid. 92, 120401 (2004) 
[6] J. Kinast et al., Phys. Rev. Lett. 92, 150402 (2004) 
[7] M. Bartenstein et al., Phys. Rev. Lett. 92, 203201 (2004) 
[8] C. Chin, M. Bartenstein A. Altmeyer, S. Riedl, S. 

Jochim, J. Hecker Denschlag, R. Grimm, Science 305 



1128 (2004). 

[9] J. Kinast et al., cond-mat/0409283 See also, Science 307, 
1296 (2005). 

[10] M. Greiner, C.A Regal, J.T. Stewart, D.S. Jin, 
cond-mat/0502411 

[11] D. Vollhardt and P. Wolfle, The Superfluid Phases 
of Helium 3, Taylor and Francis, London 1990; 
cond-mat/0012052 

[12] A. L. Fetter and J. D. Walecka, Quantum Theory of 
Many-Particle Systems, McCraw-Hill 1971. 

[13] T-L. Ho, Phys. Rev. Lett. 92, 090402 (2004) 

[14] D.M. Eagles, Phys. Rev. 186, 456 (1969). 

[15] A.J. Leggett, in Modern Trends in the Theory of Con- 
densed Matter, ed. A. Pekalski and R. Przystawa, Lect. 
Notes in Physics Vol. 115 (Springer- Verlag, 1980), p. 13. 



6 



[16] P. Nozieres and S. Schmidt-Rink, J. Low Temp. Phys. 

59, 195 (1982). 
[17] R. Haussmann, Z. Phys. B91, 291 (1993). 
[18] C. A. R. Sa de Melo, M. Randeria, and J. R. Engelbrecht, 

Phys. Rev. Lett. 71, 3202 (1993). M. Randeria, in "Bose- 

Einstein Condensation", Ed. A. Griffin, D.W. Snoke, S. 

Stringari, Cambridge Univ. Press 1995, p. 355. 
[19] J. Stajic, J. N. Milstein, Qijin Chen, M. L. Chiofalo, M. 

J. Holland, K. Levin, Phys. Rev. A 69, 063610 (2004). 
[20] A. Perali, P. Pieri, G.C. Strinati, Phys. Rev. Lett. 93, 

100404 (2004); ibid. 92, 220404 (2004). P. Pieri, L. 

Pisani, G.C. Strinati, Phys. Rev. B 70, 094508 (2004). 
[21] Y. Ohashi, A. Griffin, Phys. Rev. Lett. 89, 130402 (2002); 

cond-mat/0410220 Phys. Rev. A 72, 013601 (2005). 
[22] M. Holland et al., Phys. Rev. Lett. 87, 120406 (2001). 
[23] L.M. Jensen, |cond-mat/041243T| 

[24] J. Carlson, S-Y. Chang, V. R. Pandharipande, K. E. 

Schmidt, Phys. Rev. Letts. 91, 50401 (2003). S-Y. Chang 

et al., Phys. Rev. A 70, 043602 (2004) 
[25] G.E. Astrakharchik, J. Boronat, J. Casulleras, S. 

Giorgini, Phys. Rev. Lett. 93, 200404 (2004) 
[26] H. Heiselberg, Phys. Rev. A 63, 043606 (2001); J. Phys. 

B: At.Mol.Opt.Phys. 37, 1 (2004). 
[27] E. Timmermans et al., Phys. Lett. A 285, 2 28 (2001) 
[28] G.M. Falco and H.T.C. Stoof, cond-mat/0409621 



[29] J. Thomas, priv. comm. 

[30] G. Baym and C.J. Pethick, Landau Fermi-liquid Theory, 

Wiley & Sons, 1991. 
[31] L.D. Landau and E:M. Lifshitz, Fluid Mechanics, Second 

edition, (Pergamon, New York, 1987). 
[32] C.J. Pethick and H. Smith, Bose-Einstein Condensation 

in Dilute Gases, Cambridge Univ. Press, 2002. 
[33] H. Hu, A. Minguzzi, X-J Liu, M.P. Tosi, Phys. Rev. Lett. 

93, 190403 (2004). H. Heiselberg, Phys. Rev. Lett. 93, 

040402 (2004). 
[34] G. M. Bruun and H. Smith, cond-mat/0504734 
[35] A. Griffin and E. Zaremba, Phys. Rev. A 56, 4839 (1997). 
[36] D.S. Petrov, C. Salomon, and G.V. Shlyapnikov, 

|cond-mat/0309010| 

[37] H. Heiselberg, cond-mat/0409077 

[38] S. Cowell et al., Phys. Rev. Lett. 88, 210403 (2002). 

[39] S. Stringari, Europhys. Lett. 65, 749 (2004). 

[40] L. P. Gorkov & T. K. Melik-Barkhudarov, JETP 13, 1018 

(1961); H. Heiselberg, C. J. Pethick, H. Smith and L. 

Viverit, Phys. Rev. Letts. 85, 2418 (2000). 
[41] T. Domanski & J. Ranninger, Phys. Rev. B 70, 184513 

(2004). 

[42] D.-W. Wang, M.D. Lukin, E. Demmler, 
|cond-mat /0410494 



